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Abstract. An analytical microscopic theory for the resonant multiple scattering 
of light by cold atoms with arbitrary internal degeneracy is presented. It 
permits to calculate the average amplitude and the average intensity for one- 
photon states of the full transverse electromagnetic field in a dilute medium 
of unpolarized atoms. Special emphasis is laid upon an analysis in terms of 
irreducible representations of the rotation group. It allows to sum explicitly 
the ladder and maximally crossed diagrams, giving the average intensity in the 
Boltzmann approximation and the interference corrections responsible for weak 
localization and coherent backscattering. The exact decomposition into field 
modes shows that the atomic internal degeneracy contributes to the depolarization 
of the average intensity and suppresses the interference corrections. Static as 
well as dynamic quantities like the transport velocity, diffusion constants and 
relaxation times for all field modes and all atomic transitions are derived. 



1. Introduction 

Wave propagation in disordered media has been an active field of research for more 
than a century. The first descriptions of electronic conduction in weakly disordered 
metals by Drude M and of light propagation in interstellar clouds by Schuster Q were 
based on Boltzmann-type transport equations. Despite their successful predictions 
(e.g. Ohm's law), these theories overlooked the subtle role of interference. P. Anderson 
first showed that sufficiently strong disorder can induce a metal-insulator transition 
for quantum-mechanical wave functions on a random lattice, a phenomenon baptized 
Anderson localization S. Furthermore, it was realized that interference influences 
transport even far from the localization transition and that also other waves are 
concerned. For example, the constructive interference of counter-propagating light 
amplitudes in the backward direction, as first pointed out by Watson ||, survives 
an average over disorder, giving rise to the coherent backscattering (CBS) peak 
H [| 01- ln the past 20 years, the field of mesoscopic physics has developed, with 
beautiful experimental and theoretical results pertaining to the weak localization 
regime 

Atoms, as natural realizations of highly resonant identical point scatterers, have 
been proposed as an ideal sample for multiple light scattering and even localization 
| p| O . However, as first discovered in experiments on CBS by rubidium atoms jT^, |l4| 
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and recently confirmed on strontium atoms Jig ], the quantum internal structure of 
atomic scatterers has a strong impact on light transport properties. The Zeeman 
degeneracy of the probed dipole transition reduces the interference contrast of coherent 
backscattering dramatically. The single and double scattering contribution to the 
CBS peak for atoms with arbitrary internal degeneracy have been calculated in [ p~7| . 

In the present contribution, we develop an analytical theory of multiple resonant 
one-photon elastic scattering in an infinite medium consisting of a dilute gas of 
cold atoms at fixed classical positions, but with arbitrary internal degeneracy. Our 
results are a first step towards a microscopic generalization of the existing theories 
for the scattering of transverse vector waves by Rayleigh point scatterers |l9| 
on the one hand and of scalar waves by anisotropic scatterers [2l[] on the other. 
The paper is organized as follows. In section ||, we define the Hamiltonian for the 
coupled system light and atoms. In section |^, we calculate the average propagator 
of a one-photon Fock state inside the scattering medium with the aid of the Dyson 
equation and the photon self-energy. In section || we define the average intensity 
for the full vector field in terms of photo-detection probability. The corresponding 
Bethe-Salpeter equation is then solved by evaluating the so-called irreducible vertex 
in the Boltzmann approximation. Technically, this is achieved by a systematic 
analysis in terms of irreducible tensor operators. In section ^, we sum the so-called 
ladder series in the static limit by applying the same method of irreducible tensors 
to the intensity propagator of transverse field modes. The interference correction 
described by the maximally crossed diagrams is then easily obtained by applying 
substitution rules. In section ||, we generalize our results to the dynamic case and 
derive analytical expressions for transport velocities, diffusion constants, extinction 
lengths and relaxation times for the different field modes and all atomic transitions. 
We conclude by possible extensions of the work. 



2. Description of the system 

2.1. Hamiltonian 

In order to analyze the role of the atomic quantum internal structure, we must treat the 
atomic internal degrees of freedom quantum mechanically. For reasons of symmetry, 
the light field is treated quantum mechanically as well. Using natural units such that 
h = c = 1, the global system "matter + light" is then described by the Hamiltonian 
H = Hq + V where 

JV 

H Q = H cm + H at = ^2 a\ e a kE + ^ H at {a). (1) 

k,e_Lfc a—1 

H cm is the Hamiltonian of the free electromagnetic field. The annihilation and 
creation operators for the field modes (wave vector fc, transverse polarization e, free 
dispersion relation uik = k) obey the usual bosonic commutation relation [ake, a} k , e ,\ = 
S e ,e'Sk,k'- For brevity, a one-photon Fock state will be noted \ke) = |0) where the 
transversality e ■ k — is understood. H at is the sum of the individual internal atomic 
Hamiltonians 

Jo 

H at (a)=aj ^ \Jcm c (a))(J c m c (a)\. (2) 

m (2 — — ,J c 
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Figure 1. Resonant degenerate dipole transition, here for J = 1, J e = 2. The 
atomic resonance frequency is noted uio, the detuning of the light probe <5 = u>— luo, 
and the natural width of the excited atomic state T. The internal magnetic 
quantum numbers m are coupled to the polarization of scattered photons (full 
and dotted lines). 



Each atom a = 1,...,N is described as a degenerate closed two-level system with 
an energy separation lvq between the ground-state with total angular momentum J 
and the excited-state with total angular momentum J c (figure Q). These levels are 
resonantly coupled by the dipole interaction 

N N 

V = ^V a = -Y, D <*-E(r a ). (3) 

a— 1 a— 1 

The individual atomic dipole operator D a acts on the Hilbert space Tt g ffi Ti c of 
the internal atomic states of atom a. Dropping subscript a, we define the reduced 
dipole operator d = D/d where d = (J C \\D\\ J)/\/2J c + 1. The matrix elements of its 
spherical components d q = e q ■ d, q = — 1,0, +1 are the Clebsch-Gordan coefficients 
( J c m c \d q \ Jm) = (J e m e \ Jlmq) according to the Wigner-Eckhart theorem p^[ . In the 
dipole interaction, the electric field operator 

E(r) = i J2 £ " (^a k£ c ik r - e k a{ E c- ik r ) (4) 

k,s.Lk 

is evaluated at the center of mass r a of each scattering atom. The field oscillator 
strength £ w = (w/2eoi 3 ) 1 / 2 is given in terms of a quantization volume L 3 that will 
disappear from physically relevant expressions in the limit of an infinite medium by 
virtue of the rule L~ 3 £) fe (. . .) i-> (2tt)" 3 J d 3 A:(. . .). 

2.2. Length scales hierarchy 

The detailed description of wave propagation in a large sample of scatterers at 
quenched random positions is such a formidable task that one rather resorts to a 
statistical description, catching the generic features through configuration averaged 
quantities (average amplitude, average intensity, etc) [^| [25). In the following, 
we will consider the limit of a disordered infinite medium where resonant point 
scatterers are distributed with constant spatial density n. The absence of boundaries 
greatly simplifies the theory since translational symmetry is restored on average. The 
importance of interference depends essentially on the hierarchy between the physically 
relevant length scales in the system, here the wave length A and the mean inter-particle 
distance n -1 / 3 . We will concentrate on the low density regime where A <C n -1 / 3 . This 
inequality, alternatively rewritten as nA 3 <§; 1, justifies a semi-classical description of 
propagation along rays between consecutive scatterers that define scattering paths 
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and their associated amplitudes. It also implies that scattering paths involving 
different scatterers are uncorrelated: since their phase difference greatly exceeds 2tt 1 
the associated interference term can be expected to average to zero. Also, recurrent 
scattering sequences (visiting a given scatterer more than once) can be neglected: 
this is the independent scattering approximation (ISA) p3[ . In this dilute regime, 
the elastic mean free path (mean distance traveled between two successive scattering 
events) is given by I — l/na where a is the total cross-section for elastic scattering. 
Since we treat resonant point scatterers where a ~ A 2 , the low density condition 
implies A <C t. wave scattering is described in the far field. On the same ground, the 
equivalent relation *J~o <C £ shows that recurrent scattering is indeed negligible. And 
finally, the low density condition implies kl ^> 1 where k — 2ir/\, which defines the 
weak localization regime. 



3. Average light propagation 

The simplest quantity of interest for wave propagation in random media is the average 
amplitude, so we first compute the configuration-averaged propagator (G(lu)) . 



3.1. The one-atom scattering operator 

In the weak localization regime, a scattered wave has reached its asymptotic limit when 
undergoing a subsequent collision. Hence a central quantity in the theory is the one- 
atom transition operator T Q (z), given by the Born series T a (z) — V a + V a Go(z)V a + . . . 
in powers of the interaction V a and the free resolvent operator Gq(z) = (z — Hq)^ 1 
The scattering amplitude associated to the elastic scattering process |i) = \ke; Jm) — > 
|/) = \k'e'; Jm'} is proportional to 

(/|T a (w + iO)|t) = (Jm'\e' ■ t(u) ■ e\Jm) ^C*-*')-^ ( 5 ) 

where the notation uj + iO indicates that retarded propagators are used in the Born 
series and that the matrix elements are taken on-shell (lj = LUk = <^h' ) ■ Note that the 
incident and scattered wave vector appear only in the exponential on the right hand 
side: in the dipole approximation, the atom indeed behaves as a point scatterer, and 
the external and internal degrees of freedom are factorized. The frequency dependence 
of the internal transition operator t(uj) is 

3 172 
2npo(uj) d + ii/2 

where T = d 2 ujQ/3ireo is the natural width of the excited atomic state, S = uj — u>q 
is the detuning of the probe light from the atomic resonance and po(^>) — L 3 uj 2 /2tt 2 
is the free photon spectral density. This form of a resonant scalar t-matrix is well 
known in the context of classical point scatterers [^7| . But as the atomic scatterer has 
an internal structure, particular attention must be paid to the tensor character of the 
scattering operator. For a given atomic transition m i— » m', the Cartesian elements of 
the t-matrix are 

tij(m,m — t(u>)(Jm \didj\Jm). (7) 

This 3x3 t-matrix can be decomposed into its irreducible components with respect 
to rotations, its scalar part (or trace), its antisymmetric part and symmetric traceless 
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part: 

1 1 1 1 

Uj = -Sijtkk + -(Uj - tji) + -{tij + tji) - -Oijtkk (8) 




In the case of a point dipole (J = 0, J c = 1), only the scalar part is non-zero, 
tij = t(u>)8ij. In the more general situation of arbitrary degeneracy, we have to 
determine the influence of precisely the non-scalar parts on the average light evolution. 



3.2. Averaging procedure 

In the statistical description of light evolution inside a disordered scattering medium, 
configuration averaged quantities (•) = Tr[(-)p a t] arc obtained by tracing over the 
matter degrees of freedom. The external degrees of freedom are the uncorrelated 
classical positions {r a }i and the average is performed by spatial integration. The 
internal average is a trace over the internal density matrix. Since scattering theory 
relates asymptotically free states, this density matrix only describes the statistical 
properties of the ground state. We further assume the total density matrix p lnt = 
p,j{ct) to be a direct product of one-atom density matrices, each proportional to 
the ground state unit matrix (complete statistical mixture, no internal correlations): 

pj{a) = 27TT^ |Jm(a)>(Jm(a)| - (9) 

m 

The free evolution of the photon state is then determined by the resolvent operator 

go(z) = (G (z)) = {z-H cm y l (10) 

where the average of the free resolvent operator Gq(z) = (z — Ha)^ 1 (which is of 
course independent of the disorder) simply projects onto the atomic ground state 
chosen to have zero energy. The full average evolution in the presence of the scatterers 
is described by the average resolvent operator (G(z)) which can be expanded in the 
Born series 

(G(z)) = g (z) + go{z) (V) g (z) + g (z) (VG (z)V) g (z) + ... (11) 

Summing all repeated interactions with the same scatterer, one obtains a series in 
terms of the one-atom scattering operators T a (z) — V a + V a Go(z)V a + . . ., 

(G(z)) = g (z)+J2 9o(z) (T a (z)) g (z)+J2 9o(z) (T a (z)) g (z) (T p (z)) g (z)+. . . (12) 

The average one-atom transition operator (T a (z)) does not depend on its index 
a since all atoms are equally distributed, so that this series can be represented 
diagrammatically as 

= +N (gs +N(N- 1) ® (8 +. . . (13) 

The right hand side features the contributions of free evolution, single scattering, 
double scattering, and so forth. Starting from the following term of triple scattering 
(not shown), this average contains second- and higher-order moments of one-atom 
scattering operators, arising from recurrent scattering by the same scatterer. As usual 
in diagrammatic expansions pif , one introduces an operator representing the sum of 
irreducible contributions, the self-energy. 
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Table 1. Definition of symbols in Feynman diagrams. 



Symbol Definition 



Name 



Equation 



i-i-i—k 



90 (z) 
(G(z)) 
(Ta(z)) 
E(«) 



free photon propagator 
average photon propagator 



average one-atom transition operator (12) 



self-energy 

link between identical scatterers 



(T a (zi)T a (z2)) average single-scattering vertex 
atomic polarization vertex 



il-.jk 



X 



iljjk 



crossed atomic polarization vertex 



(11) 



3.3. Self- energy and Dyson equation 

The average one-photon Green function (G(z)) satisfies the Dyson equation 

(G(z))=g (z) + g (zMz){G{z)) (14) 

that reads in diagrammatic form 

= + •— . (15) 

By iteration, one recognizes a geometrical series which formally sums up to 

(G(z)) = [g (z)- 1 -E(z)}-\ (16) 

The Dyson equation actually defines the self-energy £(2) whose exact calculation 
remains impossible in most cases. Nevertheless, introducing the self-energy has two 
major advantages: First, any approximate expression for the self-energy will yield an 
approximate, but non-perturbative result for the average propagator (|l6|). Second, 
the perturbative expansion of the self-energy in a power series and its truncation 
are controlled by the small parameter nX 3 . The self-energy contains precisely all 
irreducible diagrams, i.e. those that cannot be separated into independent diagrams 
by cutting a single line, 



£(*) = • = N ® +N(N - 1) <g>- <g> -<E> + ... (17) 

Here, the dotted line identifies the same scattering operator appearing twice. 
Substituting this series into the Dyson equation ( fl5| ) indeed reproduces the average 
propagator ( |l3| ) in the thermodynamic limit N, L — > 00 at constant density n = N/L 3 . 

For a dilute medium n\ 3 <C 1, the independent scattering approximation (ISA) 
amounts to the first order truncation S(z) s» N(£) — N (T a (z)). With the one- 
photon matrix elements of the self-energy are 

(fcV|E(w)|fce) = N (£' • t{u) ■ e) int (e^^^^cxt. (18) 

The external average gives (e l( - fc_fc ^' r ) ox t = 5k,k>, indicating that the self-energy 
is diagonal in momentum space as required by the statistical invariance under 
translations. The internal average of the scattering operator (^) with the scalar density 
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matrix (g) is elementary using the closure relation of Clebsch-Gordan coefficients [g2j , 
(k'e'\S(w)\ke) — NMjt(ui) (e' ■ e) 8k,k>- The scalar product of polarization vectors, 
e' ■ e — S e e i shows that the self-energy, as a second-rank tensor, is proportional to 
the unit matrix, = as required by statistical invariance under rotations. 

Note that 

depends directly on the number density n = N/L 3 so that the thermodynamic limit 
is trivial. For convenience, we define the multiplicity ratio Mj = (2J e + 1)/3(2J + f ) 
with the non-degenerate limit Mq = 1. 

In optics and atomic physics, the response of an atom to an external electric field 
is often expressed in terms of the atomic polarizability ol{lj) = — (t(ui)) p8|| . In 
the ISA, the polarizability is simply proportional to the self-energy, £(<j->) = —^na{u>), 
and thus a scalar. The dilute regime is alternatively characterized by n\a\ <C 1. 

3.4- Effective medium 

Just like the free propagator and the self-energy, the average retarded photon 
propagator is diagonal in momentum and polarization: (k's'\ (G(lo)) \ks) = 
(G(k,u>)) 5 k M> (s' • e) with 

(G(k,u)) = 1 . (20) 

uj — k — h(oj) 

Its singularities are the solutions of the complex dispersion relation u> — k — = 0. 
The average propagation can be seen to proceed through an effective medium with the 
complex- valued frequency-dependent refractive index n r = k/u> = 1 — Y,(lu)/lu. The 
imaginary part of the self-energy defines the elastic mean free time 

-=-2Im£(cj) (21) 

T 

of a pure one-photon state \ke) inside the medium. Our perturbative approach proves 
to be consistent: the correction 1/r to the free frequency to is small by definition 
in the weak scattering regime 1/wr <C I. A propagating wave-packet is therefore 
exponentially damped on average with a scattering mean free path I = t (remember 
c = I). This depletion is not caused by absorption, but by elastic collisions into 
initially empty field states. At low densities, and by virtue of the optical theorem, the 
mean free path is simply £ — 1 /no~ in terms of the total elastic scattering cross section 

a = -2L 3 Im(fc £ | <I») |fee) = Mj ^ 1 + ^ 2/r2 ■ (22) 

In all of the above expressions, only a modest dependence on J and J c arises through 
the multiplicity ratio Mj. Obviously, the quantum internal structure has almost no 
impact on the average amplitude: under an average over the scalar density matrix (^|), 
only the scalar part or trace of the t-matrix (^) can survive. We are left with a scalar 
theory for the average amplitude, describable in terms of the polarizability alone. But 
one should not conclude prematurely that the internal structure has no impact on the 
average intensity which, of course, must be carefully distinguished from the square of 
the average amplitude. 
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4. Average light intensity 

In order to determine the average population of initially empty field modes, we now 
turn to the average intensity, defined in terms of photo-detection probability. 

Intensity propagation kernel 

Since the total Hamiltonian H = Hq + V is time-independent, the density matrix p 
of the coupled system "atoms + field" evolves according to p(t) — U(t)pW(t) where 
the forward time evolution operator U (t) is the Fourier transform of the retarded 
propagator G(ui) = (id — H + i0) _1 , 

1 f°° 

U(t) = -— duG(w)e- lut , t>0. (23) 

The measurable average light intensity at position r and time t is proportional to the 
average photo-detection probability , 

(I(r,t)) = AT Tr [p(t)E ( -\r) ■ (r)]. (24) 

Here, the factor Af contains the detection efficiency, and E^(r) arc the 
annihilation and creation components, respectively, of the electric field operator 
(|J). Their normal ordering assures that the photon vacuum state yields zero 
intensity. Using (p3|), we define the Fourier transform of the intensity (J(r,t)) = 
L~ 3 Eq( 27r ) -2 I duJ dn (I(q, w, 0)) e- i( - nt+q - r '> with 

(I(q,u,Q)) =M I d 3 r Tv[pG f (uj^)E^ (r)-E^ (r)G(uj + )}e iq - r .(25) 

The amplitude evolves with the retarded propagator G(v+), and the conjugate 
amplitude with the advanced propagator G^(uo-). Here, uo is the average evolution 
frequency while f2 is the frequency difference: lu± = lo ± 57/2. For small frequencies 
— > 0, the stationary regime (or long-time limit) is recovered. 

In the following, we restrict our theory to low intensity light fields, neglecting the 
non- linear response (i.e. saturation) of the atomic dipole transition, by studying the 
evolution of a field state containing at most one photon. We thus consider an initial 
density matrix of the form p(0) — p at <g> p v where p v — J2i 2 /°i,2 1 1) (2 1 describes the 
one-photon initial light field with the short-hand notation f = ki£\. For example, 
Pi, 2 = <5i,i^2.i describes a pure state consisting of an initial plane wave |fej£j). Using 
, a straightforward calculation gives the average intensity for any incident field 

(I(q, w, fl)) = N ]T pia £2 £ 3 (e 3 ■ e 2 ) 5 k3 - k2 , q $({fce}; w, SI) (26) 
1,2,3,4 

in terms of the intensity propagation kernel 

$({fc £ }; W ,0) = ((k 4 e 4 \G^(u J ^)\k 3 e 3 )(k 2 e 2 \G(u J+ )\k 1 e 1 )) (27) 
or, in operator form, $(w, f2) = (G'(uj-) ® G(w+)). 
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4-2. Bethe-Salpeter equation 

For the average amplitude, the average propagator (G{u)) had been calculated by 
solving the Dyson equation with the help of the self-energy. In close analogy, the 
intensity propagation kernel $ obeys the Bethe-Salpeter equation 

fi) = (GV_)) ® + (G f (w_)) <g> (G(w+)) I7(w, fl) fi). (28) 

Here, the irreducible vertex U(u>, Q) contains all diagrams with at least one vertical 
connection between the direct amplitude (upper line or "particle channel") and the 
conjugate amplitude (lower line or "hole channel") po| , 



U 



<g-> -*g — <g— (g) 

'>C + i + 
— <g (g> <g— X— X> 



(29) 



Factors of order iV are omitted for brevity here. One defines the reducible intensity 
vertex R(ui, Q) by 

fi) = (G t (w_))(g)(G(w + ))+(G t (w_))«)(G(w + )) fi) (G t (w_))(X(G(w + )) .(30) 

Up to a dressing by average propagators, it suffices then to solve the Bethe-Salpeter 
equation for the reducible intensity vertex, 

R(u, fi) = U(w, fi) + U(u, fi) (G f (w_)) ® (G(w+)) fl). (31) 

Diagrammatically, this equation reads 



R 




U 


+ 


u 




R 



(32) 



In general, this equation cannot be solved exactly, and the irreducible vertex U(u>, f2) 
has to be approximated. Any approximation of U(u>, f2) should be consistent with 
that of S(cj) since, physically speaking, the depletion of in initial state is caused by 
scattering into initially empty field modes. Mathematically, this consistency is assured 
order by order in perturbation theory by a Ward identity Q . 



4-3. Boltzmann approximation 

In the weak-scattering regime, the irreducible vertex may be approximated by the first 
term on the right hand side of (29) that represents the single-scattering contribution: 

U (1 \{ke};uj,n) = N ((k 4 e 4 \Tl(^)\k 3 e 3 ) (fc 2 £ 2 |r Q (co + )|fe 1 £ 1 » .(33) 

This is the so-called Boltzmann approximation (or first-order smoothing approxima- 
tion). It proves to be consistent with the independent scattering approximation ( |l8| ) 
for the self-energy since the corresponding Ward identity reduces to the optical theo- 
rem (^2|) for fl = |H . The Bethe-Salpeter equation ( |32| ) then acquires by iteration 
the familiar and simple ladder structure R(ui, Q) s» L(ui, f2), or diagramatically: 

f . . . (34) 



R 



<g— g eg— — g— _ 

Here both the direct and the conjugate amplitude are scattered by exactly the same 
scatterers. Thus, it is the squared amplitude or intensity that propagates from 
scatterer to scatterer, and all interference has disappeared. This approximation leads 
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for q — > to a Boltzmann transport equation for the intensity, justifying the Drude 
model for the electronic conductivity or the radiative transfer equation for diffusive 
light transport p3| . 

The average single-scattering intensity ([33]) can be calculated explicitly. Using 
expression (j|) for the matrix element of the one-atom scattering operator, one 
can immediately compute the external average and verify the total momentum 
conservation ( e 1 ( fe i- fc 2+fc3-fc4)-» , a) cxt — <5fe 1 +fc 3 ,fe 3 +fc 4 as required by the statistical 
invariance under translations of the infinite medium. The internal average of the 
squared atomic scattering operator ([?]) can be calculated for an arbitrary internal 
degeneracy using the techniques of irreducible tensor operators (for details, see [171] ): 

((£4 • *t(w_) • £ 3 )(£ 2 • t(u+) • £l)) int = Afjt(w_)i(w + ) I(£l, £ 2 , £ 3 , £4). (35) 

Here, I(si, £2, £3, £4) is a vertex function connecting four vectors according to 
I({x}) = wi (xi ■ x 2 )(x 3 ■ x 4 ) + w 2 {xi ■ x 3 )(x 2 ■ x 4 ) + w 3 (x 1 ■ x 4 )(x 2 ■ x 3 ) (36) 
The weights of the three possible pairwise contractions are 

so — S2 S 2 - Si Si + s 2 
wi = — - — , w 2 = — - — , w 3 = — (37) 

where the coefficients 

s K = 3(2J C + !){ ) j ^ J (38) 

are proportional to squared 6J-symbols or Wigner coefficients that are known 
to describe the possible re-coupling of four vector operators. A diagrammatic 
representation of the atomic four-point vertex has been introduced, 

l-p-,-2 1 — 2 1 2 12 

I({x})= = Wl +w 2 X + w 3 \ I- (39) 

4-LL3 4 — 3 4 X x 3 4 3 



The irreducible vertex ( |33[ ) is therefore given by 

C/ (1) ({fe£}; u, CI) = u(u, CI) _ 1- |T 2 8 kl+k3M+ki (40) 
where u(oj, CI) = N Mji(cu^)t(Lu + ) has the dimension of an energy squared. 
4-4- Weak localization corrections 

Finding among all possible diagrams the dominant interference corrections to the 
Boltzmann approximation proves to be a delicate subject | 0T| . Langer and Neal 
J32"| introduced the so-called maximally crossed diagrams yielding an interference 
correction to the ladder terms independently of the density of scatterers: 



C 



X + ">\< +■■■ (41) 



These diagrams describe amplitudes that propagate along the same scattering paths 
but in opposite directions. Their interference is responsible for the weak localization 
corrections to the conductivity of electrons in weakly disordered mesoscopic systems 
[1 and can be implemented self-consistently through U{uj, CI) w U^~'{oj, CI) + C(u>, CI) 
[30|. This interference also gives rise to the coherent backscattering peak scattered 
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from a bounded medium ||. There, the interference correction is calculated using the 
heuristic prescription R(uj, Q) ~ L(lu, Q) + C(w, Q). In sections || and||, the full ladder 
and crossed contributions will be calculated, preparing the ground for the calculation 
of the CBS peak and weak localization corrections. 

The survival of an interference term despite an average over random realizations 
may seem miraculous at first glance. Closer analysis reveals that it is really due to 
the fundamental symmetry of time reversal invariance. Let us briefly recall the exact 
relation between ladder and crossed diagrams. Take the first non-trivial term of the 
ladder sum (|34|), the second-order scattering contribution with matrix elements 

l-*^x^™^*^>-»-2 

L 2 ({k,e})= \ \ =^C/W(l,u,d,4)(G(u))(G(d))[/( 1 )(u,2,3,d), (42) 
where u = fe u e u and d = fcd£d and where the dependence on u>, fl is understood. 



Since each intensity vertex (40) conserves the total momentum, the sum also does, 
and the ladder matrix element can be written L 2 ({k,e}) = dk^ks^+k^ L 2 (q, {e}) 
with q = ki — fc 4 = k 2 — &3 and 



L 2 (q, {e}) = u{u, n) 2 }^l( £l ,e u , e d , e 4 ) <G(u)) (G(d)) l(e u , e 2 , e 3 , e d ) S kd ^ q . (43) 

u,d 

The corresponding crossed diagram 

C 2 ({k,e}) = (44) 

4-Hgi— ^!gM-3 

can be treated in the same manner, and we have to calculate 

C 2 (q c , {£}) = u(w, n) 2 ^X(£i, £ U , £3, £d) (G(u)) (5(d)) I(e Uj £ 2 , £d, £ 4 ) 4 d , qc -fe u (45) 

u.d 

where its argument qc = k\ + A;3 is now the total momentum. 

The comparison of expressions (^) and ( [45| ) shows that the value of any 
crossed diagram (the generalization to arbitrary scattering orders is evident) can be 
obtained without further calculation from the corresponding ladder contribution by 
the substitutions 

q = fei fc 4 i ^ qc = ki + fe 3 , (46a) 
(£ 3 , £4) ^ (£4, £3), (466) 

(w 2 ,w 3 ) 1 ^ (103,^2) (46c) 



The first substitution rule ( 46a ), well known in the case of scalar wave scattering, 
implies that the crossed and ladder terms are equal for scattering of plane waves 
(k\ = ki = k, k 2 = fc 3 = k 1 such that q = 0) in the backwards direction k' = —k 



since then also qc = k + k' = 0. The second substitution rule (46?), known for 
vector waves, restricts the equality to the channels of preserved helicity or parallel 
linear polarization where £"' = £. If these two conditions in the case of isotropic point 
scatterers are satisfied, then the reciprocity theorem, stemming from time-reversal 
invariance, indeed justifies the equality of ladder and crossed series. Symbolically, 
the crossed diagram can be disentangled by turning around its lower line and becomes 



topologically a ladder diagram |3J . The third substitution rule ( [46 q ) is new and arises 
because of the scatterers' internal structure: by turning around the lower amplitude 
according to the two previous rules, the atomic internal vertices (p9h are twisted. 
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To obtain the correct expression, the vertical and diagonal contractions have to be 
exchanged such that T(-, -,x 3 , x 4 ) i— > T(-, -,x±, x 3 ), or symbolically 

As explained in detail in [ p7[ , the inequality W2 ^ W3 that causes a difference of 
crossed and ladder contributions even for parallel polarizations must be attributed to 
the antisymmetric part tW of the scattering operator (represented by the coefficient 
si in equation (p7|)). 



The substitution rules ( |46a|) -( [i6q ) permit to obtain the crossed contributions 
immediately from the ladder contributions or vice versa. In principle, one can 
calculate diagrams of arbitrary scattering order using the above prescriptions. A 
diffusive transport for the intensity, however, only emerges in the limit where all 
ladder diagrams are summed up. 



5. Summation of ladder and crossed series 



Waves can be either scalar or vectorial, and point scatterers can be either isotropic 
or anisotropic. This distinction defines four classes of multiple scattering theories 
with growing complexity. The first case of scalar waves and isotropic point scatterers 
can be considered well understood |2^, ^7|. The case of scalar waves and anisotropic 
point scatterers has been solved in the framework of the radiative transfer theory p0| | . 
The multiple scattering of electromagnetic vector waves by point dipole scatterers, 
from the first approaches based on the diffusion approximation || [55], up to 
the exact solution of the radiative transfer equation by the Wiener-Hopf method 
Ji8| , pj| , has kept a somewhat discouraging appearance. Indeed, for a vector wave 
like light, polarization and direction are linked by transversality. In order to describe 
the evolution of the light intensity in three dimensions, one needs to manipulate tensors 
of rank four or 9 x 9 transfer matrices. The strategy employed in the literature consists 
in applying the scalar methods to this transfer matrix, which needs to be diagonalized 
in an appropriate way. The published results demonstrate the difficulty of the problem 
and the complexity of its solution. 

A solution for the most difficult case of vector waves and arbitrary scatterers is 
still lacking to our knowledge. But the scattering of light by atoms with a quantum 
internal structure falls precisely into this last class of difficulties since the internal 
structure couples to the polarization and the isotropic dipole approximation J = 
is forbidden by definition. In a previous work |p^| , we have been able to obtain the 
atomic intensity vertex 2({x}) for arbitrary internal degeneracy J > by a systematic 
analysis in terms of irreducible operators with respect to the rotation group. It is 
therefore natural to apply the same powerful tool in order to simplify the summation 
of the multiple scattering series as much as possible. In this section, we will consider 
the static case 51 = 0, postponing the discussion of dynamic effects to section 0. 



5.1. Strategy of summation 

Let us first recall the solution for the case of scalar waves. This case can be viewed 
as a limit of the full vector case by averaging over the incident polarization e, and 
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summing over all final directions k' and polarizations s' . The irreducible vertex ( |40| ) 
then becomes the scalar factor u = nna/ L 3 lo 2 , and the ladder series (Q) sums up to 

L(q) = tto (1 + -4(^) + A(q£f + ...) = YT^P) • ( 48 ) 

Here, the momentum transfer function A is given by the auto- convolution 

A(q£) = ^ (G(k; w)) (G(\k - q\;u)) (49) 

k 

of the average scalar propagator (^oj). The average real-space propagator (G(r;ui)) = 
L~ 3 J^k w )) c lk r in scalar form reads 

(G(r; u))=—^- e i ( fe - s ^» r (50) 
zirr 

after a suitable regularization of the UV-divergence and up to near-field terms of order 
(fcr)~ 2 [[I?]]. Using ImS(w) = —1/21 and the Plancherel-Parseval relation, one finds 

A{ql) = ¥ [ d 3 r ^ e** = ^ (51) 
v ; 47T 7 r 2 <^ v ; 

Note that no stems from the single scattering vertex f/W whereas I is defined through 
the imaginary part of the self-energy. This distinction should be kept in mind when 
approaching the localization threshold, but in the present dilute regime we can use 
no£ = 1. A thus only depends on the reduced momentum p = ql. Its small momentum 
behavior A(q) = 1 — p 2 /3 + 0(p 4 ) assures that the summed ladder propagator ( ^8|) 
has the usual diffusion pole at the origin, L(q) ~ 3uo/p 2 . 

The average intensity for vector waves is described in terms of four-point 
diagrams that connect the incident to the scattered polarization vectors. The 
vector ladder series thus defines a ladder tensor Lijki(q) according to L(q, {e}) = 
u(u>, 0) £i,i£2 j£3,fc£4,i Lyki(«jr) (here and in the following, the summation over repeated 
Cartesian indices is understood). The factor u(ui, 0) = 3na / (47r po(to)) makes the 
ladder tensor dimensionless. The ladder series ( |48| ) in tensor form reads 

L(qr) = 1 + 1 A(g) + I A 2 (q) + . . . = I (1 - A(q))- 1 (52) 

where the transfer tensor A, as in ([l9|), is the product of the atomic single 
scattering intensity vertex I and the auto convolution G of transverse propagators, 
diagrammatically 



Kjki = (Gl)ijfei = i • (53) 

: -k 



In the multiple scattering series, the tensors are multiplied in the "horizontal" direction 

(A2) ^ = ' 1 1 [ ■ ( 54 ) 

In order to stress this important feature, we group the indices by pairs left-right 
Au-jk = Aijki so that the horizontal tensor product reads explicitly 

(AB)^/;^ = Ail; mn &rnn\jk • (55) 
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The strategy of summation will be to diagonalize the tensors of rank four with respect 
to this product. We will thus try to decompose the atomic scattering vertex and the 
transfer tensor, 

i^Vr^. a = J>t<«, (56) 



in terms of suitable orthogonal projectors T^'T^ ' = JjjwT^. In this form, the 
summation of the ladder series would be trivial, 

L = ^i T W. (57) 

The corresponding crossed sum would then be obtained by substracting the single- 
scattering term and by applying the substitution rules ( f46c| )-(46c). This strategy 
proves to be successsful up to minor complications to be discussed below. 

5. 2. Irreducible eigenmodes of the atomic intensity vertex 

Let us start by analyzing the atomic ladder vertex defined by 2({e}) = 
£i,i£2,j£s,k£4,l hjkl or, equivalently, 



hi-,jk = =w 1 S i jS k i + W2 6 ik 5ji + w 3 5ii8jk. (58) 



The identity for the tensor product (p5|), a scalar with respect to rotations, can be 
decomposed into its irreducible components with respect to the pairs of indices (il) 
and (Jk) as 

i — j 

Ul;jk = = SijSki = Ys T Sh ( 59 ) 

I — k K 
where the scalar, antisymmetric and symmetric traceless basis tensors 



(60) 



define an algebra of orthogonal projectors, TWT' X ' = 8kk'^ K ^- The irreducible 
cartesian components of any rank-two tensor M are obtained by the projection 
M^- = T^ n M m „. Accordingly, we define the pairwise irreducible components 

of the atomic vertex as \( KK ) = T^IT^ \ Because the atomic vertex is globally 
invariant under rotations, only its diagonal components K = K' are non-zero, 

i = E |1 " ) = E AkT(k) - ( 61 ) 

K K 

The corresponding eigenvalues are given in terms of the contraction weights (|37 

Ao = Wl + W2 +3^3 = 1, \\ = W\ — W2, A2 = W\ + W2-(62) 



T (0) 






T (l) 


= ^(SijSki - 


- SikSji), 


T (2) 
1 il;jk 


= ^(SijSki - 


\- SikSji) — -8u8jk 
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Table 2. Eigenvalues of the atomic ladder vertex, (p3) and (pa), as a function of 
the ground state angular momentum J. 





Je = J + 1 




Je = J 


Je = J - 1 




1 




1 


1 


Ai 


J + 2 




1 


J- 1 


2(J + 1) 




2J(J + 1) 


2J 


A 2 


(J + 2)(2J + 


5) 


4J 2 + 4 J - 3 


(J-l)(2J-3) 


10(J + 1)(2J- 


M) 


10J(J + 1) 


10J(2J + 1) 



The atomic vertex function (pq ) now reads 

X({e}) = J2 x k [e^f [£ 2 £ 3 ]f >. (63) 
if 

We see that the eigenvalues e [0, 1] determine how faithfully the irreducible 
component [£164]'^ of incident field polarizations is mapped on average to [e 2 e 3 ]^ K \ 
An eigenvalue Xk — 1 means perfect mapping, an eigenvalue Xk — means total 
extinction. Note that the scalar eigenvalue is identically Ao = 1 for all J,J C . The 
scalar field mode being the intensity, this sum rule reflects the conservation of energy. 
The explicit dependence of eigenvalues on the ground-state angular momentum J is 
given in table ^| on page [l^ and displayed in figure || on page [l7| In the case of the 
isotropic point scatterer (J = 0, J c = 1) all eigenvalues saturate, Xk — 1, as expected 
for the identity operator (|5^). For J > 0, the non-scalar eigenvalues lie below unity, 
A1.2 < 1. This is consistent with the physical intuition that an initially well-polarized 
light beam will be depolarized by so-called degenerate Raman transitions between 
different Zeeman-sublevels \Jm') 7^ \ Jm). 

The "horizontal" coefficients Xk can be expressed in terms of the original 
coefficients sk defined in (|38|), 

Ak=£(-) M '(2^' + 1){ I \ K k )sk: (64) 

K' ' 

Indeed, the coefficients sk had been defined for the "vertical" coupling scheme (ij)(kl), 
and the irreducible recoupling of vector operators is a linear transformation involving 
so-called 3nj-symbols. The relation (|^) can be recognized as a particular form of the 
Biedenharn-Elliott sum rule p2| which permits to write 

X K = 3(2J C + 1) ( J ) K \ . (65) 



J C Jq J 

For the crossed series, we have to find the eigenmodes of the twisted vertex 

x^=*X J = £x* T £l (66) 

l- u -k K 

By overall invariance under rotations, the same basis tensors T^ as for the ladder 
vertex appear. The crossed eigenvalues are obtained from the ladder eigenvalues by 



applying the exchange rule (46c), 



Xo = u>\ + 3w 2 + w 3 , xi=wi-w 3 , X2 = w i+w 3 . (67) 
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Table 3. Eigenvalues of the atomic crossed vertex, ( ]6T| ) and (pq), as a function 
of the ground state angular momentum J. 
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= J + l 




Jo = J 


Jo = J 


- 1 


XO 




1 




J 2 + J - 1 
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(J+ 


1)(2J + 


1) 


J(J+1) 


J(2J4 


'1) 


XI 




1 







-1 






2J+ 1 




2J + 


1 


XI 


6J 2 


+ 12J + 


5 


2 J 2 + 2 J + 1 


6J 2 - 


1 


5(JH 


-1)(2J4 


-1) 


5J(J+ 1) 


5J(2J- 


n) 



Tabic ^ on page fl6| c ontains their explicit dependence on J, and their behaviour is 
displayed in figure |2| on page [l7|. For the pure dipole scatterer (J = 0, J c = 1), all 
ladder and crossed eigenvalues coincide trivially = \k = 1. But contrary to the 
ladder case, the crossed scalar eigenvalue xo is not fixed by any conservation law. It 
indeed deviates from unity, xo < 1) as soon as J > 0, signifying a loss of contrast 
for interference corrections to the Boltzmann intensity. In the limit J c = J — > oo, 
crossed and ladder eigenvalues take pairwise equal limits, signifying a perfect contrast 
of interference, but non-negligeable depolarization as from a classical, anisotropic 
scatterer (think of a small oriented antenna). A negative crossed eigenvalue, e.g. 
Xo = —1/3 for J c = J = 1/2, implies an even greater loss of contrast since then the 
summed series behaves like 1/(1 — xk) < 1/(1 — Ixif I)- 

Instead of using the exchange rule (1^2,^3) 1— * (11)3, W2), the crossed eigenvalues 
can be obtained from the ladder ones by a partial recoupling of vectors in the vertex. 
By using the defining properties of 6j-symbols, one finds 

X K=Y f ( 2K ' + 1 ){ I \ ?}^-S(V. + l){ 1 J Jo I (68) 

From a conceptual point of view, the expressions (|65| ) and ( |68|) of the atomic vertex 
eigenvalues in terms of 3nj-symbols are fully satisfying since we have reached their 
most concise, truly irreducible formulation. 

The direct product of two polarization vectors has nine independent components. 
The atomic scattering vertices therefore can be represented by 9 x 9 transfer matrices 
whose eigenvalues are Xk and xk- Since the atomic scattering vertices are averages 
over a scalar density matrix and therefore invariant under rotations, these eigenvalues 
are (2K + l)-times degenerate (corresponding to the Clebsch-Gordan decomposition 
of the direct product of representations of dimension 3x3 = 1 + 3 + 5). Introducing 
the irreducible components greatly simplifies a problem that at first glance defied 
an analytical treatment. It is thus natural to apply this strategy to the problem of 
transverse propagation as well. 

5.3. Irreducible eigenmodes of the transverse intensity propagator 

The sum over the polarizations of an intermediate photon in scattering diagrams like 
( |l2"|) defines a transverse projector J2e±k £ i^j = % — kikj. The average retarded 
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Figure 2. Irreducible eigenvalues of the atomic intensity vertex as a function of 
the ground-state angular momentum J. Left: ladder eigenvalues Right: 
crossed eigenvalues (pq). The unit ladder eigenvalue Ao = 1 reflects energy 
conservation. For J > 0, antisymmetric and symmetric traceless modes are not 
conserved Ai 2 < 1, reflecting depolarization. The crossed scalar eigenvalue \0 
plunges rapidly for J > 0, J e = J + 1, implying a drastic loss of interference 
contrast. In the semi-classical limit J e = J — > 00, the interference is re-established 
since xo - * 1- 
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propagator for the transverse field in momentum space therefore can be written 

(g« (*;*))= \~ fc w V (69) 

to — k — h[0J) 

Up to near-field terms, the average real-space propagator is given by (cf. (|50|)) 

(G ZJ (r;u,))=-£-S k -^A l3 (70) 

where Ay = <5y — fifj is the projector onto the plane transverse to the direction of 
propagation. Generalizing the scalar case (^) , we now need to analyze the transverse 
intensity propagator between scattering events, 



Hl;j 



fc(9)= 8^/ d3r2 ^ A - Afciei " r - (71) 

Here, the prefactor u(lo,0) = 3na / '(4tt pq(lu)) from the atomic vertex has been 
incorporated in order to manipulate a dimensionless quantity. All information 
on the vector character is contained in the direct product A^-A^ of transverse 
projectors. Contrary to the rotation-invariant atomic vertex, the intensity propagator 
now depends on the momentum q. Thus isotropic tensor modes can no longer be 
sufficient, leading to a more involved, but still exact decomposition (for details, see 



Appendix A) 



G, /:.„•,;</:• = hip) PijPki + §i(p) [PijQki + QijPki] + &s 3 (p) QijQki 

+s 3 (p) [PijPki + 2 perm.] + s 4 (p) [PijQki + 5 perm.] (72) 

in terms of the projectors Qy = cjiqj and = 5y — Qy onto the direction q = q/q 
and the plane perpendicular to it. Indeed, the projectors Qy and Py are the 
natural objects to deal with transversality and yield the most concise expressions. 
The terms in the second line of the right hand side of ( [72] ) are totally symmetric 
with respect to any permutation of indices which is indicated by "+nperm.". The 
functions s a (p) are given in equations ( A.13| ) and (A. 22) of the appendix; p — qi is 



the reduced momentum. The transverse projector in real space A^A^i has the very 
simple structure of a direct product which is obviously not conserved by the Fourier 
transformation. In momentum space, the (difficult) real-space integral equation of 
multiple scattering becomes a (simple) geometrical series of diagonal operators, but 
the price to be paid is the complicated coupling between momentum and polarization 
expressed by ( |72). 

Analyzi ng (72[) in ir reducible left-right components G^ K ' K - 1 = T^GT^ ) as 
described in Appendix A| , the transverse intensity propagator takes the form 



G(g) =£>«(P) 7f\q)+m (t (0 < 2) (<7) +T^°Hq)) • (73) 

K,a 




-0.2 1 ■ 1 ■ 1 ■ 1 ■ 1 ■ 1 

1 2 3 4 5 

P 

Figure 3. Irreducible eigenvalues jr^j ) of the transverse intensity propagator 
as functions of the reduced moment p = q£. At the isotropic limit p = th e 
eigenvalues are (2K + l)-fold degenerate, and the coupling g(p), equation (|76|), 
between scalar and symmetric traceless modes vanishes. 



The nine p-dependent eigenvalues of the 9x9 transfer matrix are partially degenerate 
such that only six irreducible eigenvalues are relevant, 

arctan(p) p 2 

~P =1 ~^" 



9 oo(p)^A(p) = ^^>=1-"- + 0(p% 



3(l-A( P )) 1 3 P 2 4 
9M= 2p* = 2-J0 +O{p) > 



3(-l + (l+p 2 )A(p)) 1 p 2 



-9(1 + p 2 ) + (9 + 12p 4 + 5p A )A(p) 7 29 P 2 4 
9M = V = 10 " 210 + ° iP } ' 



(74) 



. 6+ P 2 -3(2 + p 2 -p i )A(p) 7 13p 2 ,. 

aM = ^ = io-^r + 0(p) ' 

-3 + 7p 2 + 3(1 - P 2 ) 2 A(p) _ 7 23p 2 



The (2K + l)-degenerate eigenvalues at the origin, g#- Q (0), have been used by 
Akkermans et al. Q in a qualitative evaluation of polarization effects on the coherent 
backscattering peak, based on the diffusion approximation for scalar waves. The 
expressions at order p 2 (with a small error for the coefficient of 520) have been found 
by Stephen and Cwilich ||, (7.4)-(4.9)] and MacKintosh and John ||, (4.16)]. The 
above functions, noted \i{p) = 1 — gi<a{p), are identical to those found by Ozrin Jl8| , 
(3.16)]. In a slightly modified form, they appear also in the radiative transfer theory 
solved exactly by Amic, Luc and Nieuwenhuizen [|l9| (3. 8), (3. 15)]. But these last 
approaches do not make use of an explicit analysis in terms of irreducible components, 
which here gives also the projectors on the corresponding eigenmodes in terms of 
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Py = Sij - qiqj and Q tj = faqf. 



1 0,U;jk — 3 °ji°jfc) 

T? 3 (g)«yfc = i(PyPu - PikPji), 

T 2 1] (q)il;jk = -^{PijQkl + QijPkl - PikQjl - QikPjl), 

~To\q)a-,jk = ^(Pii - 2Qu){P jk - 2Q jk ), 

(2) 1 
T l Wihjk = -^(PijQkl + QijPkl + PikQjl + QikPjl), 

if 3 (g)«.jfc = ^(P,P« + PifcP^/) - ^P«P/*- 

For the "horizontal" tensor product (|55|), the above tensors are orthogonal projectors, 

Tf'Tr ' = b~KK'0~aifta i " 1 ■ At the limit q = 0, they recombine to give the isotropic 
tensors ( |60| ) encountered in the decomposition of the rotation-invariant atomic vertex, 
J2 a tL^"^) = T^). But this decomposition does not yield a totally diagonal 
representation. As already pointed out by Ozrin |Q, there exists an irreducible 
coupling term between the scalar and one symmetric traceless mode that vanishes 
at the origin, 



1 



(75) 



V2[3- {3+p 2 )A( P )} _ V2 



2 



m = V "" =~^f + 0(p 4 ) (76) 

with coupling tensors 

T (0 ' 2) (q)^; jfe = ^ *j (Pf* - 2Qjk), 7 (2 ' 0) (<iU:jk = ^ (Pa - 2Q a ) S jk . (77) 

The multiplication table (A. 30) of these coupling tensors can be found in the appendix. 
This coupling at q can of course be formally diagonalized, at the price of rather 
unwieldy expressions. At the limit of the summed multiple scattering series, this 
coupling term is more simply taken care of by the diagonalization of a 2 x 2 matrix. 
Furthermore, it vanishes in the diffusion approximation. 



5.4- Sum of the ladder and crossed series 

The tensor A entering into the ladder series (^) is given as the product of transverse 
intensity propagator and atomic intensity vertex, Ax,(q) = G(q) I. Knowing the 
respective decompositions (^lj) and ( |73| ) into eigenmodes, the product is simply 

Ai(q) =X>*<7* Q (p)TW +x 2 g( p )f(W +Ao5(p)f( 2 -°). (78) 

K,a 

The message of this expression is clear: in order to generalize the multiple scattering 
of light to atoms, it suffices to multiply the vector eigenfunctions gKa(p) by the 
corresponding atomic eigenvalue \k{J,J c ) — a rather trivial prescription once the 
meaning of the word "corresponding" has been clarified (which constitutes the main 
achievement of the present theory). 

The sum of all ladder diagrams is L(q) = I (1 — A.t(<jf)) -1 . To invert (1 — Aj,(<jr)), 
we try the Ansatz 

(1 - A L )-' = J2 a Ka Tf) + a 2 f (°' 2 ) + 5 f ( 2 -°). (79) 

K,a 
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By definition, (1 - A L (q))(l - A^q))- 1 = 1. The identity being 1 = Y^k,* 7 ™ > 
the coefficients of (|79| ) are completely determined. For the decoupled modes, i.e. all 
except (K = 0,2, a = 0), the result simply is a^ a = (1 — \kQko) ^ ■ The coefficients 
for the coupled modes are determined by two independent linear systems, 

£)-*-'(;)■ (£)-*-'(!)■ (8 °» 

Here, the 2x2 coupling matrix is 

a = ( 1 - , -y g ) , (si) 

\ -Xog 1 - A 2 g2o / 

with determinant |A| = (1 — Ao<?oo)(l — A2<?2o) — AoAay. The coefficients of the summed 
series are therefore obtained by inverting this matrix, yielding 

000 " |A| ' a20 - |A| ' °° ~W ° 2_ W' ( 3 
Regrouping all terms, we obtain the summed ladder series as 

L(q) = Ajf Q (p)TW (q) + A(p)T(q) (83) 

where the coefficients are 

. _ . 1 - Agggo . _ Ai _ Ai 

A oo = Ao — rn — > A n = 1 ^ > A i2 — 



|A| l-Aiffn' 1-Ai3i2 : 

A _ \ 1 — Aoffoo A _ A 2 _ A 2 

A20 = M rn ) A 2i = ; , A22 = 



(84) 



|A| 1 - A 2 .92i ' 1 - A 2 g22 

in terms of the atomic ladder eigenvalues Xk(J,J ) defined in ( p5| ) and the 
eigenfunctions gKaip) of transverse propagation (equations (|74|)). The tensors T^(q) 
carrying the information about the angular dependence have been defined in ( |T5| ) . The 
irreducible coupling term 

A(p)^A 2 A ^M ( 85 ) 



is proportional to the coupling function (|76| ) which vanishes as p — > 0. Its associated 
tensor is T(q) = ~n 2 '°)(q) + ~n°' 2 )(q) in terms of the left and right coupling tensors 
given in (|77|). 

Finally, we have to contract the external polarization vectors with the tensors, 
yielding the weights of the different field modes, 

*S W) = £ M £ 2,i £ 3 ,fc £ 4 j T<f ) (q) U ;jk > , oc . 

~ ' - _ _ ~ ( 86 ) 

^l(q) = £2,j £3,fc £4J T(q)u.jk. 

Here, the populations of the different field modes are determined by the choice of 
polarization, and are in particular independent of the atomic internal structure which 
enters only in the corresponding propagators A(p; J). The summed transverse atomic 
ladder propagator in momentum space reads 

L(q) = u(w)(j2*K a (p)t™(q)+A(p)h(4j). (87) 

K,a 

Proceeding similarly, one obtains the summed crossed tensor 
C(q)=J2XK«(p)Ti K \q)+X(p)T(q). 

K,a 
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The same tensors as for the ladder series appear, and the coefficients are obtained 
by substracting the single-scattering contribution and replacing the atomic ladder 
eigenvalues Xk by the corresponding crossed eigenvalues Xk-, 

= XO TT71 1 ) ! A ll = i > A l 



, X \ ) 1-Xi5ii l-Xi.912 , on . 

v _ , 1 - Xoffoo A v _ Xi9i\ v _ X2.922 

-^20 = X2 TTTi 1 ' A 21 = ■: , ^22 



\ X \ J 1 ~ X2ff21 1 ~ X2.922 

The determinant of the crossed coupling matrix is \X\ = (1 — Xo<?oo)(l — X2<?2o) — 
XoX2<? 2 - The coupling coefficient reads accordingly 

X(p) = X2Xo^ (90) 

The eigenvalues Xk(J, Jc) of the twisted atomic intensity vertex have been defined in 
equation (p8|) . The weights of the different modes are determined by contracting the 



tensors with the external polarization vectors (respecting the substitution rule (46?)): 
42 (<?) = ei,< £2 j £4,k £3,1 T^tfhuk, ^ 

tc(q) = £l,i £2,3 £4,fc £3,i ~t(q)iijk- 

As for the ladder case, the population of a field mode is determined solely by the 
choice of polarization vectors (but is different from the corresponding ladder mode 
unless £4 = £3). The atomic internal structure only enters in the corresponding 
propagators X(p; J). The total summed crossed propagator therefore reads 

C{q) = u{Lu)(^2x Ka { P ) tg2W) + X(p) t c (q)) . (92) 

6. Transport velocity, diffusion constant and relaxation times 

Let us now turn to the discussion of dynamic quantities like the transport velocity 
and the diffusion constant for the different field modes. In the dynamic setting 
!] ^ 0, the atomic ladder and crossed vertices ( |6l| ) and (^6|) remain unchanged. 
Indeed, the only dependence on frequency had been factorized in the prefactor 
u(u), f2) = N Mjt(uj + )t(uj-) and included in the transverse ladder kernel (|7l|), 



Gu- Jk (q,n) =u(tu,n)L 3 J d?r{G lJ {r-u J+ ))(G kl {r-Lo^))e l "- r (93) 

where u>± = lu ± 0/2. Then, 

G ihjk (q,Q) = J d r — — — A M e iq r (94) 

where the frequency dependence of the propagators gives rise to an ^-dependent 
complex scattering mean free path 

m = -w , h( \ -o - (95) 

1Zj(CJ + ) — 1Zj(U>_) — 111 

The static limit is of course £(0) = I = -l/2Im£(w). Using E(w) = NMjt(w), the 
prefactor of the integral (j94|) can be written 



B{n) - iehp — r (96) 
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such that 23(0) = 1. Obviously, the tensor structure and q-dependence in ( |94| ) are 
identical to the static case (]7l|), so that the Fourier-Laplace transforms yield the same 
results, now featuring the complex mean free path £{Q). The eigenfunctions of the 
transverse propagator become 

g Ka (q,n) = B(Q)g Ka (p(Q)), (97) 

where the static eigenfunctions gKa(p), given in (|74|), are evaluated at the dynamic 
reduced moment p{£l) — q£(Q). Since the atomic eigenvalues Xk and \k are not 



affected, the ladder and crossed sums are still given by (87) and (|92|), with the 
dynamic eigenfunctions ( j97[) and the dynamic prefactor u(w,tl). Please note that 
these expressions are exact in fi and q. 

6.1. Transport velocity 

In order to determine the long-time and long-distance behaviour of transport, one 
habitually develops the propagator to lowest orders in Q and q. Since the atomic 
internal structure has completely factorized from the frequency dependence, we 
exactly recover results that are well-known from the case of resonant scattering of 
scalar waves by point particles [^3[ Let us briefly show how the main results 

are obtained very easily in our framework. The small frequency behavior of (p6|), 
0(f2) = 1 + iOrtr + 0(il 2 ), defines a common transport time scale for all field modes, 

E' (uj) 

r tr = Im-ff +*(1-Re£»). (98) 

where E' = dE/dw. Its two contributions have simple physical interpretations. The 
first term on the right-hand side can be traced back to the prefactor u(u>, fi) and is 
simply the (Wigner) time delay 

Tw = Im xR = ^ (99) 

where 4>(u>) is the phase of the scattering t-matrix t(u>) — Yj(uo)/NMj. The second 
term is the propagation time £/v g between consecutive scatterering events with the 
group velocity inside the effective medium, v g = du/dk = (1 — Re E'(w)) -1 (remember 
c = 1 and the real part of the dispersion relation, k(uj) — u> — ReE(w)). The 
group velocity by itself looses its physical meaning in the vicinity of a scattering 
resonance where extinction cannot be neglected. The transport velocity w tr = £/r tI = 
(WgT 1 +rw) _1 , however, stays causal, v tr < 1, at all frequencies u), and may decrease by 
orders of magnitude at resonance for a high enough density of scatterers [ 39 , 23| . In the 



present case of resonant point scatterers, the above expressions simplify remarkably 
since 

ReE'M E'M 1 , . 

+Im—^4- = (100) 

2ImE(w) E(w) r V ' 

up to terms of order (S/oj) 2 , (T/lu) 2 . The transport velocity then takes the simple form 
«tr = (1 + (r^) -1 ) -1 , which corresponds to an inverted Lorentzian of width L because 
of the reso nant character of the inverse scattering mean free path € _1 oc l/(5 2 + T 2 /4). 



The sum (100) has been called "dwell time" (in the case of non-resonant (L — > oo) 
scattering of electrons, the dwell time vanishes due to a Ward identity that is invalid for 
scattering by frequency-dependent potentials |39|] ). The total time of transport ( p8|) 
can be written as the sum of the free propagation time £ between consecutive scatterers 
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and the dwell time 1 /Y. At resonance, for example, the Wigner time delay is t\v = 2/T, 
but the group velocity correction is negative since Re£'(cj )/2ImI](ajo) = — 
Whereas the Wigner time delay accounts for the (possibly large time) spent "inside" 
a resonant scattering object, the total dwell time includes the self-consistent dressing 
of scatterers by the surrounding effective medium which must not be neglected. 



6.2. Extinction lengths 

The transversality of the light field and the atomic internal structure are connected to 
the spatial variable q. Setting fl — in a first step, all propagation eigenfunctions ( |r4| ) 
have a quadratic development in p = q£ of the form gKa(p) ~ ok~ CKaP 2 , whereas the 
coupling term vanishes quadratically, g(p) = cp 2 . At order p 2 , the coupling between 



scalar and symmetric traceless modes disappears, and the propagators (84) for the 
transverse ladder modes behave like 

Here appears an extinction length 



e Ka (J) = tJ (102) 

V A K — "K 

that is responsible for an exponential decay of the mode K, a in real-space (in field 
theory, this corresponds to a finite particle mass m = £ Ka )- We notice immediately 
that the intensity or scalar mode K = has a diverging extinction length since 
bo = 1 and Aq = 1 for arbitrary internal degeneracy. This true diffusion pole (or 
massless Goldstone mode) survives thanks to a fundamental conservation law, the 
conservation of energy. All other field modes have finite, and in fact rather short 
extinction lengths £ko.{J) ^ A implying that the transversality of propagation mixes 
well-defined polarization modes. The smaller \>k and Xk, the shorter the extinction 
length. Thus, the atomic internal degeneracy, responsible for A1.2 < f as soon as J > 0, 
leads to even smaller extinction lengths for non-scalar field modes. This observation 
confirms the intuitive picture that degenerate Raman transitions between different 
Zeeman sublevels | Jm!) ^ \Jm) contribute to scramble the field polarization. 

Analogous arguments apply for the crossed propagator. Here, the atomic ladder 
eigenvalues are replaced by their crossed counterparts \k, and the crossed extinction 
lengths are 



Sk«(J)=IJ T\ ■ (103) 

The scalar eigenvalue \o 1S n °t constrained by energy conservation and can become 
very small for J > (cf. figure ||). The atomic internal degeneracy breaks the 
time-reversal symmetry between ladder and crossed propagators, and leads to a rapid 
exponential damping of the interference modes. 

Extinctions lengths of the order of or smaller than the scattering mean free path 
£ imply that the evolution can no longer be considered diffusive. Amic et al. jl9| 
have stressed that the exact extinction length, defined as the inverse of the pole of 
the propagator with smallest imaginary real part, can differ from the above diffusive 
expression by as much as 50%. This difference is bound to become even worse as the 
atomic eigenvalues decreases. Clearly, a quantitative prediction for the propagation 
of the full vector field must go beyond the diffusion approximation. This is especially 
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true if one is interested in coherent backscattering from a finite scattering medium, 
where short paths or non-scalar field modes can become dominant. 



6.3. Diffusion constant and relaxation times 

Combining the developments linear in O and quadratic in q, the ladder propagators 
take the form 

h-Ka (<Zi Q) « [ b ^l 1 - ■ (104) 

—l\l + t K + Djiaq 

Here appears the diffusion constant D^a — (ScKa/bfc) Dq for the corresponding field 
mode, simply proportional to the diffusion constant of the intensity, Dq = ^i> tr /3, 
and independent of the atomic eigenvalues. Furthermore, a relaxation time has been 
defined, 

tK = n /v-i i (105) 

The physical meaning of this relaxation time is equivalent to that of the extinction 
lengths £fc a defined above since l 2 Ka = D^a^K- For the scalar mode of field intensity, 
t^ 1 = implies a truly diffusive transport as required by local energy conservation. 
For the non-scalar field modes, the finite relaxation times t± 2 become of the order 
of the transport time t% t indicating that light transport can no longer be described 
accurately in the diffusion approximation, a tendency aggravated by the internal 
atomic degeneracy leading to Ai,2 < 1. For the crossed propagator, analogous 
conclusions hold. The quantum internal structure introduces dephasing times as will 
be discussed elsewhere |Bo| . 



6.4- Conclusion and things to be done 

In summary, we develop a consistent theory for the multiple scattering of photons 
by a dilute gas of cold atoms. The external degrees of freedom of the atomic point 
scatterers are supposed to be classical Poissonian variables. Particular attention is 
paid to the internal degrees of freedom: a resonant dipole transition of arbitrary 
degeneracy is treated analytically by a systematic use of irreducible tensor operators. 
We sum the ladder diagrams of the full transverse vector field, and calculate the 
correction of maximally crossed diagrams. The internal degeneracy has no impact on 
the properties of the average amplitude (such as the scattering mean free path) since 
the average over a scalar internal density matrix projects onto the scalar component 
of the transition matrix. Furthermore, the internal degeneracy is only coupled to the 
(static) polarization vectors and completely factorized from the (dynamic) frequency 
dependence. Therefore, the transport velocity and the diffusion constant remain 
unaffected. However, the interference properties are strongly modified. Indeed, the 
non-scalar parts of the scattering t-matrix survive in the average intensity and are 
responsible for a decrease of contrast as soon as J > 0. In the diffusion approximation, 
we give expressions for extinction lengths or relaxation times that depend in very 
simple manner on J. Atoms thus appear as an important class of anisotropic scatterers 
with intriguing interference properties that allow a complete analytical description of 
light scattering. 

The present contribution deals with transport inside an infinite scattering 
medium, taking full advantage of the statistical invariance under translations which 
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makes all operators diagonal in momentum representation. The influence of boundary 
conditions, relevant for coherent backscattering or transmission experiments, may be 
evaluated in the framework of the exact Wiener-Hopf method or of the approximate 
method of images. Both approaches deserve a detailed discussion which is beyond 
the scope of this paper and will be discussed elsewhere pi] ]. Another important 
issue is the influence of an external magnetic field. Perhaps the systematic use of 
irreducible tensors can help to simplify the theoretical description of light transport 
in magneto-active media . In addition, the impact of a splitting of the internal 
Zeeman degeneracy remains to be studied, experimentally as well as theoretically. 

The transversality of the propagating light field imposes a change of polarization 
in the course of multiple scattering. In this respect, the present theory provides a 
microscopic analogue of spin-orbit coupling studied extensively in electronic disordered 
systems [|| . But contrary to the electron case of spin \ , no weak anti- localization for 
the spin 1 photon can be expected, however strong the spin-orbit coupling. The 
scattering of photons by atoms with internal degeneracy appears as an analogue of 
spin- flip scattering of electrons by magnetic impurities Q . Interestingly, we can derive 
exact expressions for characteristic extinction lengths and relaxation times, and further 
experimental as well as theoretical studies can be envisaged. The links between optics 
and atomic physics on the one hand and condensed matter physics on the other thus 
promise to continue to be an important source of inspiration to both fields. 
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Appendix A. How to find the irreducible eigenmodes of the transverse 
intensity propagator 

Appendix A.l. Decomposition in real space 

Starting from the purely transverse projector in real space Ga-jk{v) = A^Afc;, where 
Ay = Sij — rifj, let us first determine its "left" and "right" irreducible components 
by projecting onto the isotropic basis tensors (|6C|), 

G C*.*0(#) = TWG(r)T(*'>. (A.l) 
The exchange symmetry (i, j) *-* (I, k) implies that the sum of orders K + K' is even 
since G^^f \f) = Q>\fj,f \f) = (-1) K+K ' G^-* } (r), by parity properties of the 
j{K) _ This condition decouples the antisymmetric from the symmetric modes. We 
find the purely scalar component 

= I Wik, (A.2) 
the antisymmetric component 

G$$(r) = ±{A ij A u -A ik A jl ), (A.3) 

the traceless symmetric component 

1 12 

G ii-jk( f ) = o( A io A ki + \k^ji) - ~(5 u A jk + 5 jk Au) + - S u 5 jk , (A A) 
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as well as two scalar-symmetric mixed components 

The irreducible components of the transverse propagator G(q) are thus given by 
integrating the previous expressions according to (|7l|), 




(A.6) 



We introduce the short-hand notation G^ KK '(q) = f q & K K '(f) for this Fourier- 
Laplace transform. The angular integral of the expressions ( |A.2| )-( |A~5| ) depends on 
the number of times that the components fi of the unit vector appear: we distinguish 
scalar terms of type 5ij5u, quadratic terms of type fifjSki and a quaternary term 

fifjf k fi. 

Appendix A. 2. Scalar component 

The angular integral of the scalar terms is trivial, of course, and the result of the 
Fourier-Laplace transform is 

jl = s (j>) = \A(qt), (A.7) 

proportional to the scalar transfer function A(p) — arctan(p) /p. The scalar component 
of the transverse propagator therefore is simply given by 

G^^iq) =A(q£)T^ (A.8) 

using the definition ( ]60j ) of the scalar basis tensor T^ ). Naturally, the scalar 
component only depends on the transfer function already known from the purely 
scalar case. 



Appendix A. 3. Antisymmetric components 

The angular integral of the quadratic terms is easily solved using a generating function 
argument, 

?^e- = -4-5-^M. (A.9) 
47r r z oqidqj qr 

Call for brevity s(x) = sin(x)/a;, x = x = rq, di = d/dxi, and s'(x) = (d/dx)s(x) 
so that (A.9) evaluates to 

- d% s(x) = -s" (x) Mj - ^ (Sa - Mj ) • (A. io) 

The derivation (or equivalent angular integration) generates the quadratic terms q\q^ 
that could be expected, but also a new term proportional to Sij. Define the two 
orthogonal projectors on the subspaces parallel and perpendicular to q, Qij = q^j 
and Pij = Sij — qiqj. They add up to the identity Sij and satisfy 

^im^mj ^ijl QimQmj Qij ) PivaQraj QimQmj 0- (ATI) 

The Fourier-Laplace integral of quadratic terms is thus given by 

fifj = §t(p) P i:j + s 2 (p) (A. 12) 
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where the coefficients are 



^..if^.^)- 3[1- (1+^)] 



2p J x 4p 2 

3 r^„-x/.^_3[i-^(p)] 



O /.CX3 

S 2 (p) = -± dxe-*^s"(x) = 
2 PJo 



(A.13) 



2p 2 



The three functions hip) (defined in (A. 7)), §i(p) and hip) are not independent. 
Indeed, the co ntrac tion of indices in the quadratic term ( A. 12] ) must reproduce the 
scalar integral (A/7), so that §o(p) = 2s\(p) + h(p)- 

The purely antisymmetric component of the transverse propagator (|A.3| ) contains 
only scalar and quadratic terms, so that the previous results permit to write 

G (M) (<?) = 9ii{p)T^\q)+g 12 (p) T^(q). (A.14) 
The projectors depend on the direction q, 

T{ 1} (q) tl ;jk = \(P l3 Pkl - P l kP J l), 

- , (A ' 15) 

(4)il;jk = 7:(PijQkl + QijPkl ~ PikQjl — QikPjl)- 



Thanks to the relations (A.ll), they are indeed orthogonal projectors for the 
horizontal tensor product (|55|), T^tL 1 ^ = S a pTa . By construction, they are purely 



1 P 

antisymmetric from left and right, so that the antisymmetriser T*- 1 ) commutes with 
them, T^Ta^ = Ta^ — T^T^. The eigenvalues of the decomposition ( A.14 ) depend 
on the reduced momentum p — q£, 

3[1 - A( P )] 



gulp) = hip) 



912 (?) 



2p 2 

_ s (p) - h(p) _ 3[-l + (1 +p 2 )A(p)} 



(A.16) 



2 Ap 2 

It is instructive to consider the limit of zero momentum q — > 0, where all dependence 
on the direction q must vanish. Both cigenfunctions 

1 ^tj 2 1 n 2 

9u(p) = 2 - !J + °(P% 9i2(p) = g - + 0(p% (A.17) 

have the common limit fjia (0) = \ , so that the two q-dependent tensors recombine to 
the isotropic antisymmetric projector, 

T^(q)u ;j k + Tf\q)*vk = ~ hkh) = T 4 ( g fc . ( A - 1 *) 

The antisymmetric mode K = 1 must have three eigenvalues that are degenerate 
at zero momentum where indeed gi a (0) = h- At finite momentum p > 0, the 
degeneracy is partially lifted, and the two distinct eigenvalues gn(p) and g\2 (p) 
appear. Ozrin obtains the same eigenfunctions (noted 1 — Xg(p) and 1 — As(p), 
respectively, (3.16)]) and shows that the remaining twofold degeneracy is 

carried by the function 312 (p)- Note that we express the propagator in Cartesian 
components G^'^f ' and not in its re-coupled spherical components G^^ ' = 

Y^ pqrs {Km\llpq){K'm'\llrs)G^rs ' which turn out not to diagonalize the propagator 
and therefore are not well adapted for our purposes. 
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Appendix A. 4- Symmetric traceless components 



For the symmetric traceless components ( A. 4 ) of the transverse propagator the angular 
integral over a quaternary term, 

47T 

has to be calculated, giving 



■r i r j r k r l e i «- r = dt jkl s(x) (A.19) 



1 / s \ , „ „ „ , / s 



dfjkMx) = - - (P li P H +2perm.)+ - (P«Qw+5 perm.)+s""(x) Q t] Q kl . (A.20) 

The result must be totally symmetric with respect to any permutation of indices 
which is indicated by "+nperm.". The total Fourier-Laplace transform then takes 
the following form, 



fifjf k fi = s 3 (p) (P i:j P k i +2 perm.) + s 4 (p) (PyQ fci + 5perm.) + s 5 (p) QijQ kh (A.21) 
with the p-dependent coefficients 

2pJ x \x J 16p 4 



S4 (p) , A f d.e-^ (fX = 3 + 2^-3(1+^) (A.22) 
2pJ \xj 4p 4 



s 5 ( P) =ir^-* /p s""(*)=- 3+p2+3A{p) 



2pJo ' ' W 



Again, these functions are not independent because a contraction of indices in (A.21) 
must reduce to ( A.12 ), implying §i(p) — 4s 3 (p) + s 4 (p) and s 2 (p) — 2s 4 (p) + s 5 (p). 

Knowing (A. 14), the symmetric traceless components of the transverse propagator 
can be predicted to have the form 

G (2 ' 2) (S) =92o(p)T^(q)+ 921 (p)j[ 2 \q)+g 22 (p)T { 2 2 \q). (A.23) 

Indeed, for zero momentum q = 0, the three anisotropic tensors T„ (q) must 
recombine to the isotropic projector T^- fe = ^(SijS k i + Si k Sji) — \o~uo~jk- At non-zero 
momentum, the degeneracy is lifted, and each identity is replaced by S rs = P rs + Q rs , 
such that 

T (2) _ 2 n n 

' il;jk ~ g Wil^jk 

+ \{PnQki + QijPki + P lk Q 3 i + QikPji) - \{PuQj k + QaPjk) (A.24) 
+ \(PijPki + PikPji) - \PaPjk 

By applying the projector T^ 2 ) to the terms of order 0, 1, 2 in separately, we find 
the following traceless symmetric expressions, 

T ( o } (q)iiuk = \{2Qa - Pu)(2Q jk - P jk ), 

(q)a;jk = \(PijQkl + QijPki + PikQjl + Q lk P 3 i), (A.25) 
1 { 2 2) {q) a , jk = -(PyPu + PikPji) - \PaPjk- 
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(2) (2) (2) 

Again, these tensors are orthogonal projectors, T„ TJ,' = 8 a p T„ ' . By construction, 

(2) (2) 

their partial left and right traces vanish, [Ta) m m;jk = {Ja )u-,nn = 0. The 
corresponding eigenfunctions are 

2 -r^-^i9-^ -9(l+p 2 ) + (9 + 12 P 4 + 5p 4 M(p) 



,,_ 3 3 6 + p 2 - 3(2 +p 2 -p 4 ).4(p) 

-3 + 7p 2 + 3(1 - p 2 ) 2 .4(p) 



322 0) = s 2 (p) + 2s 3 (p) = 

confirming the expressions for A2 = 1 — 3201 -^4,6 — 1 ~ 321 and A3. 8 = 1 — 322 of Ozrin 
pl[ eq. (3.16)]. We recall their behavior close to the origin, 

7 29p 2 ^. 4x 

920{p) = Tq-^w +o{p) > 

92 1 (p) = ^- 1 -^ + 0( P % (A.27) 
Appendix A. 5. Mixed symmetric modes 

Finally, the mixed scalar-symmetric modes of the transverse propagator are 

G (0 < 2) (S) = 3(p)T (0 ' 2) (g), G^(q) = g(p)f^°Hq), (A.28) 
in terms of the mixed tensors 

f ^ 2) {q)u; 3 k = ^ S u (P jk - 2Q jk ), l {2ft) {q)a, 3k = ^ (Pa ~ 2Q«) S jk . (A.29) 
These tensors are no longer projectors; their multiplication table is 

j(k,k>) t (k») = s K , K „8 a oT< K < K \ (A.30) 

The coupling function is given by 

g(p) = ^[3q(p)-35 2 (p)] _ y/2[-3+(3+p 2 ).4(p)] 
6 4p 2 
and vanishes quadratically at the origin, 



2 



.9(p) = -^ + 0(p 4 ). (A.32) 
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